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ABSTRACT 


In this study, steady state cooling of a mixture of CO 
and CO 2 using laser radiation has been investigated for 
industrial applications. The basic cooling scheme involves 
exciting CO molecule with a frequency doubled CO 2 laser, 
resonance energy transfer to CO 2 upper laser level followed 
by fluoresence emission to the ground state of CO 2 . The 
system configuration consists of a tube through which the 
mixture of CO and CO 2 flow while a circular laser beam 
falls perpendicular to the tube. Appropriate equations 
have been developed to calculate cooling rate in the radiated 
as well as post radiated zones. In the analysis, the effect 
of radiation trapping has also been included. 

Our results show that cooling is produced m both 
radiated and post radiated zones when the total pressure is 
less than 200 mtorr. Further, the cooling rate is a function 
of pressure and tube radius j cooling rates are lower for 
smaller tube radii. The maximum volumetric cooling rates 
are in the range of 6-90 W/M and are obtained for pressure 
90-100 mtorr and tube radii 1-2 cm. 
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CHAPTER 1 


INTRODUCTION AND LITERATURE SURVEY 

1.1 INTRODUCTION 

Low pressure gases could be cooled using a conventional 
heat exchanger wherein a coolant is passed over a container 
which encloses the gas to be cooled. In a conventional heat 
exchanger we rely on convective as well as conductive modes 
of heat transfer. For sufficiently large velocities, the con- 
ductive heat transfer is small as compared to the convective 
heat transfer. The various parameters which influence the 
convective heat transfer rate m a conventional system are the 
geometry, the Reynolds number and the Prandtl number. However, 
for a fully developed laminar flow m a channel, the heat 
transfer rate is independent of Reynolds and Prandtl numbers 
[l] . Essentially, the conventional systems are based upon 
surface cooling. 

In the recent past, several schemes have been proposed for 
cooling low pressure gases using laser radiation [2-8]. A 
majority of them propose to reduce the thermal velocity of 
atoms or molecules, thereby reducing the Doppler width. These 
schemes, therefore, have applications in high resolution spectro^ 
scopy. A different scheme of laser induced cooling has also 
been proposed in ref. [8] and this scheme has applications m 
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Photon engine. In this work, a laser cooling scheme proposed 
recently [7] has been extended for bulk cooling for industrial 
applications . 

The scheme analysed consists of pumping an excited 
state of an atom or molecule followed by near resonant 
energy transfer to another atom or molecule which then 
reradiates at higher or lower frequencies. Since the 
resonant energy exchange process can take place over the 
entire volume of the gas radiated by the laser beam, bulk 
cooling can take place. Bulk cooling has an advantage in 
achieving uniformity m heat exchange as compared to a 
conventional system where temperature differences within the 
volume are bound to occur. 

Another advantage of the scheme analysed in the present 
work is that cooling can be achieved remotely without another 
fluid being brought m contact. This may find application m 
cooling m hazardous environment as m nuclear reactors. It 
could also be applied in outer space where a laser beam housed 
in a satellite can effect cooling on another satellite or a 
planet. At some places where local cooling is required, this 
system will prove to be very effective, whereas it will be 
difficult to achieve the desired results using conventional 
systems . 



An important feature of this cooling system is that it 
is independent of the flow geometry and thermal properties 
such as thermal conductivity. Therefore, it could find 
applications m uniform cooling of a mixture of gases, m 
which the constituent gases have wide ranging thermal proper- 
ties . 

In this work, we have considered an actual system which 
will use laser radiation to cool a mixture of flowing gases. 

In this system, the laser beam is located perpendicular to 
the flow direction. The cooling effect m the radiated zone 
and the after effects have been calculated taking radiative 
trapping into account. Comparison with the conventional 
system has been made and discrepancies are discussed. Further, 
those regions over which laser cooling will be more effective 
than the cooling in conventional systems have been explored. 

In Chapter 2 the basic principle used In effecting laser 
cooling is discussed and appropriate rate equations are derived. 
In addition an expression for cooling rate which takes into 
account radiative trapping has also been developed. In 
Chapter 3, the variables on which the cooling rate is dependent 
are Identified. These variables are then optimized to get 
maximum cooling rate. In Chapter 4, the post-radiation zone has 
been explored and the heat removal in it' Is computed. Chapter 5 
contains the laser power calculations. It also shows the 
comparison of this system with a conventional system. Summary, 
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conclusions and recommendations for future work are included 
in Chapter 6. 

1.2 LITERATURE SURVEY 

In the literature one finds essentially two types of 
cooling schemes using lasers s transient cooling and steady 
state cooling schemes. The transient cooling schemes are 
based upon one of the effects : radiative recoil [2,3], 
collisionally aided Flourescence [4] and the storage of 
energy m a vibrationaly stable mode [8]. Hansch and 
Schawlow [2] were the first to propose such a scheme in which 
it has been shown that magnesium vapours at a temperature of 
600°, when illuminated by intense light on the low frequency 
side of the Doppler broadened singlet resonance line at 
2852. 1°A, would result in a reduction m the average velocity 
by a factor of 50 and a ,24°K drop m the temperature. A 

system for generating an optical trap with a capacity of 

7 a 

3x10 atoms and trapping rate ^ 10 atoms/sec, using beams 

tuned 50 y^ (natural line width) below the sodium D resonance 

frequency with a CW power 200 mV/, has been proposed by 

Ashkins [3]. This could result m average velocity ^ 3 cm/sec 

_AO 

which is equivalent to 10 K. 

The steady state cooling methods make use of multilevel 
systems and are based on : absorption of laser radiation follow- 
ed by reradiation at higher frequencies [5,6] and resonant 



energy transfer by an excited atom or molecule to another 

atom or molecule which then radiates at higher fiequencies 

[7]. In all these schemes radiative trapping has been 

ignored. Djeu [5] has proposed to use laser radiation to 

resonantly excite P branch transition on higher J number 

in a molecular gas. The excited gas flouresces on P,Q 

and R branch transitions on low J number , releasing some 

of the photons of higher energy. The difference m energy 

is taken from the translational degree of freedom. For a 

Xe rich mixture using a 1 watt laser, he has evaluated the 

temperature difference between the wall of the tube and the 

axis as &^.l r, if r is expressed in cm. Goela and 

Thareja [6] have proposed a scheme of using laser radiation 

to resonantly excite molecules having high Frank Condon 

factors but lower frequencies and getting reradiation on 

high Frank Condon factor but higher frequencies. The 

scheme is such that it could also be used for effecting 

cooling m high pressure gases. For a CN molecule, pumped 

with a tunable laser in the wavelength range 4159.6 A 0 to 

4208.2 A 0 , they have evaluated the heat removal rate as 
-13 —1 

10 cal S per excited molecule for a laser power 
3 —3 

'•"'‘10 W and pulse energy 10~ J. Goela and Thareja [7] 
have proposed another scheme of using resonant energy 
transfer to produce cooling m a mixture of gases. For a 
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mixture of CO and CO 2 at .1 torr, they have calculated a 
heat removal rate of .286 miA? using a frequency doubled CO 2 
laser of .5 mW. This scheme has also been used m this 
study for producing bulk cooling m a gaseous mixture of CO 
and CO 2 (as discussed in details m Chapter 2). 



CHAPTER 2 


ENERGY REMOVAL IN THE RADIATED ZONE 

2.1 INTRODUCTION 

In this chapter, we have outlined the procedure which 
we have followed m making energy removal calculations in a 
laser irradiated zone. First we have discussed how cooling 
effect is produced with a laser. Then a schematic represen- 
tation of the proposed system is shown. A small note on rate 
constants is followed by the derivation of the governing 
equations which will be used m later chapters for detailed 
calculations . 

2.2 LASER INDUCED COOLING 

Laser Is a source of intense energy. Therefore, one 
would expect to produce heating with it. However, under 
certain conditions, the resonance interaction between laser 
and atoms/molecules may also produce bulk cooling of the gas. 
This could be understood from Fig. 2.1. Here G and H are 
two species of atoms or molecules and 1,2,3 and 4 could 
either be electronic energy levels, vibrational levels or 
rotational levels. A laser of appropriate frequency excites 
atoms/molecules of G from level 1 to level 2. The excited 
atoms/molecules of G collide with those of H in energy level 4, 
and resonant energy transfer takes place as given below. 




Fig. Scheme of producing steady e;ct* M'p-ing 
through energy transfer 
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G(2) + H(4) — G ( 1 ) + H( 3 ) - <^E (2.1) 

where E is the energy diffeience given by 

E = (E 3 -E 4 ) - (E 2 -E 1 ) = (E 3 -E 2 ) (2.2) 

where E (l = 1 to 4) is the energy of the ith level. Both 
E ^ and E^, being the energies m the ground states, may be 
taken to be zero relatively. 

Once the atoms/molecules of H are excited to level 3, 
they decay to level 4 via spontaneous emission. For cooling 
it is necessary that (E^-E^) > (E 2 ~E^), so that A.E is 
positive m eqns. (2.1) and (2.2). Cooling of the mixture is 
observed because a part or whole of this spontaneous radiation 
leaves the system. The cooling effect is brought about by a 
decrease m tho translational energy of the gaseous mixture. 

In this work the cooling scheme for industrial applica- 
tions is centred around the above scheme. The specif will 
G will represent the gas CO and H the gas C0 2 . The level 2 of 
CO corresponds to 2143 cm~\ and level 3 of C0 2 is at 2349 cm 
Accordingly 

AE l2 = E 2 -E 1 = 4.26 xlO” 20 J (2.2a) 

and 

AE 34 = E 3 -E 4 = 4.67xl0~ 20 J (2.2b) 

As a result, 

A E 32 = E 3 -E 2 = 0.41xl0" 20 J (2,2c) 
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Since E32 > 0 , cooling 
conditions are favourable 
later m this chapter as 


can be achieved, provided other 
,as we will explore m detail 
v/ell os m Chapter 3 . 


2.3 SCHEMATIC REPRESENTATION 5 CP LASER COCLING SYSTEJ 

A flow system, m which the gaseous mixture to be 
cooled is brought to interact with the laser radiation, will 
provide steady state coolinc . The geometry as shown m Fig. 
2.2 explains the scheme further. It essentially consists of 
a tube of radius R through which the gaseous mixture 
(CO + CO2) flowo at a flow speed u. The laser beam is also 
circular m cross section havmo a radius R and it falls 
perpendicular to the direction of flow. It is obvious that 
if the laser beam has the same radius as the tube, the 
interaction between the gaseous mixture and the laser would 
be maximum. A beam with radius greater than R, would 
remain partially unutilized. 

Before we proceed to evaluate the rate of change of 
translational energy, it wn 11 be illuminating to first 
present a small note on vibrational relaxation rates and 
resonant energy transfer. 

2 . 4 VIBRAT IONAL REL \XAT ION 


Considerable experimental information Is available on 
the rates of exchange of energy be tween the vibrational 
and translational modes of a molecules. In a mixture of 



Gaseous “ 
mix turc “ 

(flow _ 
speed u)- 

(co + co 2 )- 



taser radiation 
( Circular section) 


Fig 2*2 Schematic representation of laser 
cooling system 
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two gases (Fig. 2.1), the relaxation of the excited molecules 
of G and H can be representeu by the following equations : 


G(2) 

) \/ -—£**r>* 

“T ,/Y ™ 

G ( 1 ) 

+ X 

^ E l2 

(2 

.3) 

H( 3) 


H ( 4 ) 

+ X + 

A b 34 

(2 

.4) 

where X 

is any 

nolec 

ule of 

the mixture, 

and AE l2 and A 

E 34 


are as given m eqn. (2.2a) and (2.2b). 

This, m reality, is a simplified way of visualizing 
the relaxation process which allows it to be characterized 
by two rate constants only, viz., kp^ (s -1 ) for reaction 
(2.3) and k ^ (s“^) for reaction (2 .4). Thus is m tune 
with the experimental information available on the rate 
constants, espec^'lly of the nu .fcure of CC and C02? as will 
be dealt with m the present work. 

Since X m eqns. (2.3) and (2.4) is any molecule of 
the mixture, the rate constants k 0l and k 0<1 will be considered 
as functions of the total pressure of the mixture. The 
actual reaction rates of reactions (2.3) and (2.4) will also 
depend on the partial pressures or number densities of 
G(2) and H(3), respectively. These will appear explicitly 

V 

m Section 2.6, when the rate equations are written. At 30GK, 
the values of and k^ are given by [9,10] 

I<2f 10“ 2 P T (s" 1 ) 


(2.5) 
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and 

1< 34 = 193 p T (s" 1 ) (2.6) 

where is the total pressure of the mixture m torrs. 

The rate constants k and k ^ can easily be expressed 
as functions of the total number density n^ using the ideal 
gas relationship 

p = nkT (2.7) 


where 

p = pressure m P a 
n = number density, m“ 2 

O 

k = Boltzmann constant, 1.38x10“ J/K 
T = temperature, K, 


and the conversion factor between torr and (1 torr v-%133 P„). 

3 . a 

The resulting values are : 


k 2l 


3.105xl0“ 25 n T 


( 2 . 8 ) 


and 


k 34 = 5.99xlO“ 21 n T 


(2.9) 


-1 -3 

where and k^ 4 have units of s and rij in m 


It may be pointed out here that these values- should be considered 
applicable in the vicinity of 300K. Variation with temp, will 
not be accounted for in this work. 



2.5 RESONANT ENERGY TRANSFER 


The resonant energy transfer between energy level 2 of 
G and energy level 3 of H can be represented as follows : 

G( 2 ) + H( 4 ) ■ > • G ( 1 ) + H ( 3 ) ~ A t (2.10) 

G(l) + H( 3 ) — ^G(2) + H(4) - & E (2.1l) 

where &. E is as m eqn. (2.2). 

The rate constant associated with (2.10) will be denoted 
by k 23 ( s ~^) a nd that associated with (2.11) by k ^2 (s”^)« 
These rate constants will be considered as functions of 
partial pressures of H(4) and 0(1), respectively. The 
dependence of the actual reaction rates of (2.10) and (2.11) 
on G(2) and H(3), respectively, will be written explicitly 
m the rate equations m Section 2,6. This is m tune with 
the experimental data available, where, m fact these rate 
constants are given as function of total pressures. However, 
if significant change m level population is involved, these 
must be considered as depending on the respective ground level 
populations , 

With this interpretation, the values of k 23 and k 32 are 

r 9] : 



k 23 = l* 3 8xl° 3 Pi (s" 1 ) 

(2.12) 

and 

k 32 = 3.68xl0 3 p 4 (s” 1 ) 

(2.13) 
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where and p 4 are the partial pressures in torr of molecules 
m energy levels 1 and 4, respectively. 


As 

ctions 


m Section 2.4, ^03 
of number densities n 


and k 
^ and 


32 

n 4 


can be expressed as fun- 
of the respective levels : 


k 23 = 4.285xl0~ 20 n x 


(2.14) 


and 

-20 

1< 32 = 11.42x10 n 4 (2.15) 

where k2 3 and k 32 are m units of s~ E and n^ and n^ are in 


2.6 RATE EQUATIONS 

As pointed out earlier, the cooling effect is a result 
of decrease m the translational energy of the gaseous mixture. 
The rate of change of translational energy of the gas per unit 
volume can be written as 


dt 


n 2 k 21 ( E 2" E 1^ + n 3 lc 34^ E 3" £ 4' ) 
" n 2 k 23^ E 3”' E 2^ + n 3 k 32^ E 3" E 2^ » 


(2.16) 


where 

O 

E-j. = translational energy of the gas mixture, j/m , 
n i = number density of molecules in level i, i = 1 to 4,rrf 3 , 
and k ±i are the rate constants as explained earlier. 
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The population n of various energy levels depends upon 
the laser pumping and other energy transfer mechanisms. These 
populations, for the radiated zone, may be calculated from the 
following eguations : 

dn 0 

dt“' = W 12 n l " U 21 n 2 “ n 2 lc 21 “ n 2 k ?3 + n 3 k 32 ” n 2 A 21 ( 2 * 17 ) 


dn, 

dfiT 


"‘ f '34 n 3 + l<: 23 n 2 


n l +n 2 = n G 


n 3 +n 4 = n H 


K 32 n 3 - n 


V 

3 lX 34 


(2.13) 

(2.19) 

( 2 . 20 ) 


where 


W 12 



h..i 


t 


34 


n 


G 


n 


H 


are the rates of stimulated absorption and 
emission respectively, between levels 1 and 2, s 
is tne rate of spontaneous emission from level 2 
to 1, which for CO for the levels in Fig. 2,1, 
is 33.8 s" 1 [10] 

is the rate of spontaneous emission from level 3 
to level 4, which for CO^, for the levels in 
Fig. 2.1, is 400 s" 1 [10], 

is the number density of gas G in m™° which is 
CO in the present work 

-3 

is the number density of gas H m m , which is 


CO 2 in the present work 



n H + n G 


1 f 


llso. 



( 2 . 21 ) 


where 

is the total number density (m ) as explained 
m Section 2.4. 

The eqns. (2.17) and (2.18), however, should be 'notified 
to take into account radiation trapping. This is explained m 
the next section where appropriate modifications are mtrocuced. 

2.7 RADIATION TRAPPING 

In radiation trapping, a photon emitted by an excited 
atom/molecule is absorbed by another atom/ molecule, reemitted, 
reabsorbed and so on. This reduces the effective rate of 
spontaneous emission of radiation from to FAjj [ll]. 

According to HolsLein r s analysis for cylindrical geometry and 
Doppler broadening, 


i- 


1 .60 

k Q R(it log k Q R) 


" 1/2 


( 2 . 22 ) 


where R is the tube radius (m), 

k Q is the absorption coefficient, m~\ at the center of 
the line given by the expression : 



g(o) (n ± 



(2.23) 


flic re 


.2 

rS 0 2 q » 

Bur) -1 


(2.23a) 


A. 4 


S 1 

g(o) 


is the wavelength corresponding to the transition 
is the spontaneous transition probability, 
is the refractive index of the medium, 
is the degeneracy of level i, and 
is the line shape function. 


3 


1 y 


For the case of Doppler broadening, the lane shape 
function, g(o) is given by 


2 /log 2s 1/2 


(2.24) 


Ay xs the width of the Doppler oroadenec- line (full width 
at half maximum) given as 


(AV) d = — / ( m f 


(2.25) 


where 

n is the mass of the atom or molecule, Q is the 
frequency corresponding to the transition j i, 
k is the Boltzman’s constant, 

T is the temperature of the gas m Kelvin and 
c is the velocity of light m vacuum. 

In an actual situation factors producing both types of 
broadening i.e., Doppler as well as collision broadening may be 



present at the same lime, however, at low pressures, it is 
the Doppler broadening that dominates. As will be seen 
shortly, pressures of interest to us ere sufficiently low. 
Hence, only Doopler broadening effect has been taken into 
account m the calculation of th j radiation escape factor F. 

The ramje for the values of F is from 0 to 1. The 
extreme case of F being equal to zero is when no radiation 
leaves the system. Physically, this will be the case when 
the number densities are large and/or the Lube radius is large. 
The other extreme of F being equal to one occurs when the 
number densities are very low and/or the cuoe radius is small 
and all the radiation manages to escape from the system. 

Hence, from now on, we will work with F instead of 

A . , F . accounting for the trapping of radiation correspond- 
ing to transition between the states j and 1 . 

2.8 GOVERNING EQUATIONS 

The modified rate equations can be written as 


dn^ 

dt 


vV 12 n l 


W 21 n 2 


n 2 k 21 


n 2 k 23 + n 3 k 32 


21 A 21 n 2 (2.26) 


d n 3 

dt” = ”^34^34 n 3 + Ic 23 n 2 “ k 32 n 3 “ n 3 k 34 


(2.27) 


For the steady state condi'cion, there is no change m the 
population of level 2 and 3, i.e,, 



and 


( 2 . 20 ) 


dno 
d t" 


O, 


dn. 


-t •- 

O u 


0 


Therefore, on solving the aoovc equations, we obtain 


n. 


k 23 n 2 

k 34 + k 32 + F 34 A 34 


(2.29) 


Hr 


U 12 n G 


' u if 

u 23 32 


C ( ' , 12 +¥i 21 +F 21 A 21 +k 21 +k 23 )- 


(2.30) 


where use of eqn. (2.19) has also been made. 

To explore the conditions under which cooling can take 
place, we make use of eqn. (2.29) in eqn. (2.16) to write the 
latter as : 


d t" 


n 3 (E 2 “H 1 ) . [ (F 34 A 34 -i-k 32 +k 34 ) 


k 21 , „ . E 3~ E 2 1 
k 23 + K 34" r 34 A 34 tbj- eT J 


(2.31) 


For a mixture of CO and C0 2 , as m our case, the first term in 
eqn. (2.31) involving k 23 is likely to be negligible m compa- 
rison with the other two terms. This follows from an order of 
magnitude comparison of the values given m eqn. (2.8), (2.9), 
(2.14) and (2.15). This will be further confirmed in Chapter 3, 
where it will be shown that this term is indeed negligible for 
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the C0-CQ 2 system under all operating conditions of interest. 

V/e will, therefore, make simplifications by neglectino 
this terra m all our calculations. Under these conditions, 


“ »3<V E d 


E 0 ~E 0 

f k -F A - -?.] 

L 34 ^34 a 34 E 2 -E ] _ J 


(2.32) 


dE t 

Since is the rate of change of translational energy, cool- 

dE . 

mg will be observed if < 0, i*e„, 


F 34 A 34 


p 

' i E 2 “E 1 ; ? k 34 


(2.33) 


The upper limit of the pressure at which this condition will be 

satisfied can be obtained by assuming the radiation escape 

factor F.^ = 1 and using the values given earlier, together 

with eqns. (2.6) or (2.9). The .limiting total pressure is 

found to be 0.199 torr. In terms of the total number density 

21 “3 

of the mixture, it corresponds to n^ = 6.42x10 m 

Obviously, the actual operating pressures have to be below 
this limiting value. These, together with corresponding partial 
pressures, will be explored in detail m Chapter 3. 



CHAPTER 3 


OPTLlUIu COOLING CONDITIONS 

3.1 INTRODUCTION 

In this section we will ideality the factors on which 
the cooling rate depends. The varjatxon m cooling rate as a 
function of these variables will be discussed. This will be 
followed by a brief note on the technique used to optimise 
cooling rate. Results using this technique are then presented 
in tabular as well as graphical form and subsequently, the 
total cooling effect and the resulting temperature drops are 
evaluated and discussed. 

3.2 RADIATION ESCAPE FACTOR CALCULATIONS 

The expressions for radiation escape factor F have been 
listed in Section 2.7. The relevant data and the numbers 
arrived at by using these expressions are presented below. 

For CO at 300°K corresponding to levels 1 and 2, we have, 

A 21 = 33,8 sec" 1 

* l,507xl0 8 Hz 
~ 5143 crft = 4*66 x 10~ m 
= 1.824X10” 19 [nj^] nT 1 


m 


D 


* 

k 


12 

CO 


(3.1) 
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and for CO2 at 300 °K, corresponding to levels 3 and 4, 

A3 4 = 400 sec” 1 

= 1.3177x10 s Hz 
X 

^34 =23X9 Clfl = 4.26x10 m 
CO 

lc 2 = 1.353xlO“ 18 (n 4 ~n 3 ) pT 1 (3.2) 

The values derived here are used to evaluate the radia- 
tion escape factors namely F 2 ^ and ^34* ft nay be indicated 
here tha t f 21 is not needed for the radiation zone calculations, 
it is required later m Chapter 4 for the calculations in the 
post-radia tion zone. 

3.3 IDENTIFICATION OF VARIABLES 

For the system presented m Fig. 2.2, it is easily seen 
that the cooling rate will depend on the following variables 

a) Total pressure p^ (or number density n^) , 

b) partial pressure of one of the constituent gases, 

P C0 or p co 2 ' 

c) flow speed u, 

d) tube radius R, and 

e) laser intensity I. 
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However, if we make the assumption that the characteristic 
residence time of a molecule m the radiation zone is 
sufficiently larger than the average relaxation time, the 
energy removal rate does not depend on tie flow speed, u. 

As will be seen later, a result of this is that the temp, 
drop m the raoia tion zone is inversely prooortional to u. 

The cleoendence of the cooling rale on the tube radius R 
is only through the dependence of the radiation escape factor 
F 34 on R, as can be seen from eqi. (2,22), In order to account 
for the dependence of the cooling rate on laser intensity, we 
introduce a parameter s which is a measure of saturation, as 
follows : 



The maximum value of s is 0,5, which will occur if the laser 
beam is sufficiently intense. Otherwise s < 0.5. 


Similarly, to account for the partial pressures of the 
two gases in the mixture, we introduce a parameter r defined 
as follows : 


r 



(3.4) 


obviously, o < r < 


1 
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In earns oC 'Weso p:> ca \"0 :.n /rite fro a err. 

(2.32) the follow nc c^p'^csion for 'be volume Wic cooling 
3 

rate, (w/n ) ; 


C R 


dE j 
dt 


C 23 • s * r • ( ^“31^/34^^^ 

^ F 34 A 34 + ^ 1 “ S ^ C 32 rn T +sC 23 rn f' i ' C 3e n T^ 


(3.5) 


where 034,023 and CUo are the numeric a l factors aooaaring in 
eqns, (2.9), (2.14) and (2.15), resoec Lively . It is seen th~t 
C^ depends on s,r and n~. explicitly as well as implicitly 
through F^- , while dependence on R is only implicitly through 
F 0 )/ j_. we will first discuss Ihe dependence of Cp on the para- 
meter s, which is a measure of saturation between levels 1 and 
2 as indicated above. 


3W EFFECT OF II '2 SmTURATIU i PARA. lETEK ON IKE COOLING RATE 

If F-p ,ff ero not to depend on the saturation parameter s, 
it is easily seen from eqn . (3.5) that the cooling rate Cp will 
be a mono ton ic'-lly increasing function of s. Thus the maximum 
value of Cp, w: 11 occur for the maximum value of s (i.e. 
s = 0.5), for an> value of other parameters like, r, n^ and A, 
But m reality, s affects n^ and n^ , ’which m turn affects F 34 . 
Since the algebraic dependence is too complicated, quite cum- 
bersome expressions are obtained if one pursues the matter 
analytically. No will therefore be satisfied by investigating 



.he dependence of on s, numerically. 


For this purpose, we calculate C,-, for some given values 
Df R, n.p and r, while continuously vorung s from 0 to 0.5. 

The calculation procedure is shown m the form of a flow chart 
m Fig. 3.1. i.'lost of the calculation involved is for evaluat- 
ing Fg^ iteratively. The cooling rate is finally obtained using 
eqn. (3.5.) , 

o 

The values of the cooling rate (w/m°) are plotted in 
Fig. 3.2, for R = 0.5 cm, n^ = 3x10^ n“° (pj = 93 mtorr), and 
for two different values of r. It is seen that the maximum 
cooling rates are indeed achieved for s = 0.5, i.e., at the 
saturation condition. do will mate a further comment on it 
m Section 3.7. mos t of the remain.) ng calculations m t h is 
work will be cone for the saturation condition, for which Cp 
can be obtained jy substituting s = 0.5 m eqn. (3.5) : 



2 

C23 r ( 1***^ ) n-j f F “"Cmg n-p 6H 2 i] 
[ 2F34A34K *^ 23 + ^32 ^ ■^^■j , t'2C0^np ] 


(3.6) 


3.5 OPTIMAL PRESSURE RATIO FOR A SPECIFIED TUDE RADIUS AND 
TOTAL PRESSURE 

Figure 3.3 plots the cooling rate as a function of the 


parameter r, which is the ratio of the partial pressure of gas 
G (i.e., CO) to the total pressure of the mixture. This is 
done for a fixed tube radius, ana two fixed values of the total 



Ana Numerical constants 


n^rny 


ny=( l-r)ny 

s =0.025 


4 34 


Evaluate n^ using eqn.(2.29) 
n 4^H“ n 3 

Evaluate k 2 using eqn.(3.2) 


Calculate F^* using eqn.(2.22) 


Yesj 

"1 TS ^ 
.s<0.5 


S — S“f" S 





Calculate C R using eqn.(3.5)^ 


Fig, 3.1 Flow chart for evaluating the cooling rate 

for given values of parameters,R,ny and r with 
a consistent calculation for F 34 











Total pressure = 93 m torr 



Fig32 Volumetric cooling rate at different conditions of saturation 
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prvscurn. Flo calculations are, of course, clone for the 

^nturam on condition, i.e. s = 0.5. 

It is seen that for given values of tube radius and total 
pressuie, there is an optiraal value of r which leads to maxi- 
mum cooling rates. This is also brought out from eqn . (3.5) 
or eqn. (3.o), which shows that = 0 for r = o, and also 
for r = 1, indicating one or more extremal points between 
r = o and 1. 

It should also be noted that for a given value of n T 
(or Pj) > there may be cooling for all values of r, or there 
(aav bo heating for some values of r, and cooling for other 
values of r, as is clear from Fig. 3.3. In Fig. 3.4 we also 
show the corresponding radiation escape factors. I ^ is seen 
that F 34 increases as r increases, which is to oe expected, 
because increase m r signifies decrease in 302 

We have plotted Figs. 3.3 and 3.4 in detail for illustra- 
tion only. Otherwise, our interest lies m identifying, for 
given values of R and n^, that value ox r v\fo ich will provide 
maximum cooling rate. Towards this end, we use an optimiza- 
tion scheme based on Fibonacci method [12]. This is discussed 
in Section 3.6, and the results obtained are discussed in 
Section 3.7. 
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3.6 FIBCK.CCI ,.1511- OD 

This method con be used to fnc che unium or maximum of 
t func ci on it one variable, even if t'ie funct_on is not conti- 
nuous. This method has the fo llowing limitations. 

(i) The initial interval m which the optimum lies, has to be 
known. In our case this is known because o < r < 1. 

(ii) The exact optimum can not oe located. Only an interval, 
known as the final interval of uncertainty will be known. The 
final interval of uncertainty can be made as small as desired 
by making more computations. 

This method makes use of the sequence of Fibonacci 
numbers, F n , defined as 

F = F, =1 
o i 

and Fj = F I-JL + Fj„ 2 , I = 2,3, ..., N (3.7) 

yielding a sequence 1,1,2,3,5,8,13 .... 

It essentially consists of comparing the values of the function 
at two points, that are specified using Fibonacci numbers, m 
the interval of uncertainty. A portion of this interval is 
discarded based on this comparison and the value of the fun- 
ction at a new point (again specified using Fibonacci numbers) 
m the remaining interval is evaluated. A comparison is made 
again and the process is continued till a desired interval of 



unc^r Lciini <.y is achieved. 
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The optimal can lie an/where within 
-' K : s final inta'' r al. By cl.oosing U = 20, we can calculate the 
value of r wi Ji on accuracy of better than + 0.0001. 

The flow chare is given m Fig. 3.4. Fo^e and 

are tne extreme limits of the variable r, i.e., = 0 and 

Bj_ = 1, and F = 20. The program based on this chart will 
calculate that value o> r r, i.e*, ( Aq/ n^) , for which coolinc 
rate is maximum, for a given value of tube radius (R) and 
total number density (n^). As can be seen from the chart, R 
and r\j are used as input values. 

The cooling rate Cp is given by eqn. (3.6) if s = 0.5, 
or cqn . (3.5) if s < 0.5. The scheme used to obtain a consistent 
value of Fn/ for evaluating the cooling rates is same as in 
Fig. 3.1. It nay also be pointed out that various other values 
of interest, e.j., n 3 ,n 4 etc. can also be obtained using the 
flow chart shown m Fig. 3.2. As will be seen later, these 
values are required as input conditions to the post radiation 
zone. 

This exercise has been repeated for various values of R 
and n T . The results thus obtained are presented and discussed 
in the next section. 

3.7 RESULTS AMD DISCUSSION 

Using the Fibonacci method we can obtain the maximum 
cooling rates corresponding to the various values of n-j. for 
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CR1 =C r ( r x ) 

CR2=C d ( r 0 ) 


Compare CRl ’ 'ith CR2|- 
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Fig. 3.5 Flow chart for calculating optimal value of 
r by Fibonacci method 
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A sot of values obtained using Fibonacci method 


Tube radius .5 cm 


_A-q" 

n T xlO u 
(No./m 3 ) 

P T 

(mtorrs) 

r 

(n G /n T ) 

Cooling rate 
(W/m 3 ) 

5 

15.5 

0.4884 

4.78 

10 

31.0 

0.6856 

13.10 

15 

46.6 

0.8026 

18.90 

20 

62.0 

0.8472 

22.90 

25 

77.0 

0.8738 

25.10 

30 

93.0 

0.8919 

25.80 

35 

108.0 

0.9051 

25.06 

40 

124.0 

0.9150 

23.10 

45 

139.0 

0.9229 

20.00 

50 

155.0 

0.9293 

16.00 

55 

170.0 

0.9347 

11.00 

60 

186.0 

0.9392 

5.30 

65 

202.0 

0.9990 

-0.0018 
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different tube radii, and 'he corresponding values of the 
parameter r. Table 3.1 shows the results of this exercise 
foi a tuoe ramus O.o cm. Similar values have been obtained 

for other tube radii. 

From Taale 3.1, we can easily see that maximum achievable 
cooling rale for this tube radius is achieved for n T = 30xl0 20 / 
m (or p T = 93 mtorrs) and at r = 0.892. Since at these 
values maximum cooling is achieved, we call these conditions 
as optimum operating conditions. The cooling rate attained 
with these values is 25.8 w/rrm The total pressure p-j- is 
calculated using ideal gas relationship p = nkT, for 
T = 300 K, where k is the Boltzmann constant. In Fig. 3.6, 
the values of the maximum cooling rates are shown for the 
saturated and unsaturated conditions. Once again, it is 
seen that the maximum cooling rates are achieved for the 
saturation condition. 

Figure 3.7 shows maximum cooling rate as a function of 
total pressure for tube radii 0.1 cm and 1 cm. Fig. 3.8 shows 
the same for tube radii 1 cm and 2 cm. From these two plots, 
it is easily seen that maximum achieveable cooling rate per 
unit volume is the highest for tube radius 0,1 cm. This is so 
because for low radius of the tube, the radiation trapping is 
low and a significant amount of radiation leaves the system. 




pressure for 
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Fig. 3-7 Comparison of maximum cooling rates for R-0*1 and 1cm 




=1cm 
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Fig.3*8 Comparison of maximum cooling rates for R =1 and 2cm 
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Figure 3.9 snov/s the optimal conditions p-j. and r with 
tube radius as a parameter. The optimal r increases with 
increasin'! tube radius. This is so because as the tube radius 
increases , ihe radiative traopmg increases . One way to 
decrease this effect is by reducing the amount of CO 2 gas 
present m the system. Increase m r has the effect of 
lowering the CO2 content of the mixture. This however is not 
necessarily required for low tube radii operating at low 
pressures, as is seen from the plot for R = 0.1 cm at lower 
values of p^ m Fig. 3.9. 

From representative detailed data presented above, the 
following points are clearly brought out : 

(a) Corresponding to each tube radius, there is an optimum 
total pressure of the mixture, which leads to maximum cooling 
rate per unit volume. 


(b) In order to achieve this maximum cooling rate, the partial 
pressures of the two gases should have a definite value. 


These results are given in Table 3.2. Also listed in 
this table are the total cooling rates, E R1 , obtained as 
follows : 


E 


R1 



D 3 C R 


(3.8) 


3 

where C D is the volumetric cooling rate, W/m . These 
K 

results are also shown in Figs* 3 #10 and 3*11* 
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k 21 


Oh THE COOLING RATES 


As v.vi r. pom ten out m Chapter 2, Section 2.8, the terra 

mv living k 0 ^ in the cooling rate expression has been neglected 

in ,i 1 1 our calculations (see eqns . (2.31) and (2.32) ) . It is 

easy to soo the effect of this term on the cooling rate and 

the outimal opera ting conditions listed in Table 3.2 . For 

example , for R = 1 cm, it can be seen that this term reduces 

3 

cool mg rates by approximately 0.05 W/ m , This is negligible 

3 

co ape; rod to the actual cooling rate obtained, viz., 13.5 W/m . 

Ellect of this term is of the same order of magnitude for 
other operating conditions. Hence, neglecting it is entirely 

jus ti : red. 

3.9 rEuPERATUiE DROP IN THE RADIATION ZONE 


Using the optimum cooling rates m the radiated zone, we 
will evaluate the temperature drop m this zone for various 
tube radii using the following expression : 


C(/Cp)cO + (J>C p ) c02 3 . A.u. A© 1 = E Ri C3.9J 

2 

where, A is the area of cross section of the tube (m ), 

u is the flow speed (m/s), 

is the temperature drop m the radiation zone, 
and C p are the density (kg/m 3 ) and the specific heat 
(j/kgK). The subscript denotes the species of which 
the properties must be used. 




(auun)OA nun jad) 6ui)ooo ajqDAaitjOD uunuijxD^ 


Radius (cm) 

Fig 3-10 Maximum volumetric and total cooling rates for different tube radius 
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‘jc isily c » n o ith species ( 1 ) is evaluated using the ideal 

-:us Invv : 


i 



(3.10) 


k as ‘.nil man* s constant = 1.38xl0 -23 j/°K, 

T 

R, is U.o gas constant ] and n (n~ 0 ) is the number 

1 kg°K 1 

of the ith species . 

The guS constant (R) for CO is 296.9 -™r~ 

kg K 


and for C0 o this value is 188.9 - — 

kg°K 


fhe specific heat (C ) for CO is 1113 ™ — 

p kg°K 


and for COr> it is 837 — — [ 13] . 

kg K 

Using the values of listed m Table 3.2, the values 
of u A0j_ can be calculated. Table 3.3 shows the maximum values 
of ( u A9j_) achievable in the radiated zone for some tube radii. 
It can be easily seen that a flow speed of 0.1 m/s gives rise to 
a temperature drop which is 6.9 K for 0.1 cm tube radius and 
13.5°K for 2 cm tube radius. The temperature drops for the 
remaining tube radii lie in between these two values. If larger 
temp, drops are desired, lower values of the flow speed can be 

used* 
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Table 3.3 


v;uuis of u achievable m the radiated 

a f ' a function of tube radius 


S.Vu 

Tube radius 

( Ciil) 

Maximum value of 

uAS, [fJS] 

1. 

0.1 

0.69 

n 

* m f> 

O.b 

1.17 

3. 

1.0 

1.29 

4. 

l.b 

1.32 

5 . 

2.0 

1.35 



F03T RADIATION ZONE 


4.1 3 ' ,1 Rt J.UTION 

in the I ollowmq sections, we discuss the cooling (or 
heating ) cl feet m the region succeeding the radiation zone, 
rau. *t.l snows the two zones namely the radiation zone and the 
pest radiation zone . In the analysis we will assume the 

'allowing : 

(3) The number densities or partial pressures of molecules 
m different energy levels is uniform throughout the radiation 
rene. Tnis assumption is justified because the laser beam is 
sufficiently intense leading to saturation in the entire region 
thereby giving rise to uniform number densities of molecules 
m level 1 and 2 . 

(L) The spatiaJ diffusion of molecules of a particular 
energy state, due to gradients m x-direction in the number 
densities m the post radiation zone is negligible. 

By virtue of assumption (a), we immediately know the inlet 
conditions to the post-radiation zone. Assumption (b) simpli- 
fies the balance equation to be written for the number densities 
in the elemental volume Adx shown m Fig. 4.1. 

It may also be pointed out here that one does not know, 
a priori, whether there will be a temperature drop or rise in 
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t h ■ ■ Cion zone. this is because whereas any trapping 

:,i r : r'l. i.> '■'ucly o. -lLtod radiation, as well as the collisional 

represented oy k 2 ^ and will lead to 


r* 1 iX.it i '.r. 


' O L S * 


r.i- . t . n t t ■ *' r« sonant energy transfer to level 3 and subsequent 
usee, e : i odui ion leads to cooling. Thus the overall effect 
can only hi uunu by a formal analysis cf the post-radiation 
ran *, wbach differs from that for the radiation zone in 

fall t,v j r, * re meets : 


(a) Ir the post-radiation zone there is no pumping of 

me Lucid from level 1 to 2. 

(b) The number densities are functions of x. 

A/' RATE EQUATIONS 

Vie bruin by defining the reference x = o, at the end 
0 i the red i utj.cn zone where the number densities n ± » i = 1 to 
4 ate known from calculations m Chapter 3. Although the end 
of this f onc will be a curved surface, we will assume it to be 
a piano passing through the outer most point of the radiation 
zone. One may also consider x to be measured from the curved 
exit surface of the radiation zone. 

Writing a balance equation for n 2 , assuming that there is 
no accumulation of molecules m level 2 in the volume Adx 
i.e. , assuming steady state, we get % 

Un 2 (x) = un 2 (x+dx) + (F 21 A 21 n 2 +k 23 n 2 - k 32 n 3 )dx 


(4.1) 
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,t M 1 


* i: 


* 


i 


Vm v " " " 4 7 ^ tshu? to be uniform over the 

* * * r -'’ ^ 'low area A cancels out. Also k 0 , 

z.±. 

''" c earlier. Si'tilaily for level 3, we 


Ui 3 (x) i.n ;i ( :-r-i •:) + ( r 3 4 A 34 n 3 +k 32 n 3 - k 23 n 2 + k 34 n 4 )dx 

(4.2) 

Iqr.c. ( 4 . 1 ) ar.u (4.2) can be written m the f orm of following 

di: f*'» < r-t . i] c iuations : 


t\ x 


‘‘ ::) A ?l n 2 “ k 23 n 2 + k 32 n 3 


ana 


u 


” F 34 A 34 n 3 + k 23 n 2 " k 32 n 3 “ k 34 n 4 


(4.3) 

(4.4) 


I r. troth:. mg < now variable *X, = , Hqn. (4.3) and (4.4) can 


in* y/rd kN-n 


dm, 

‘cT'C “ F 2J A 2l n 2 “ k 23 n 2 + k 32 n 3 


(4.3a) 


an., 

d'Z " F 34 A 34 n 3 + k 23 n 2 " k 32 n 3 " k 34 n 4 


(4.4a) 


On cube ti Luting proper expressions for k 1 j , we get 


dn 2 

= ~ F 2l A 2l n 2 ~ c 23 n 4 n 2 + c 32 n 3 n l 


(4,5) 


dn q 

= “ F 34 A 34 n 3 + c 23 n 4 n 2 ” c 32 n l n 3 


c 34 n 4 n T 


(4.6) 
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- V ■* * ^ a 1 \ .t % ) are two coupled differential equations 

1 c ' - ' ; using the standard Rung e-Cutta-G ill 

* ’ [ l 4 ’ 1 . , once the initial conditions are known. 

4. ’ It V. A. : . ■ ! to [TICKS 

In % t*v ' , we nave derived the operating conditions 
r. u;.vl , t. , «;/• :.-j r or various tube radii. These can easily 
bt- .1 t a . : • • i i iM.i Table 3.2. In eqns. (4.5) and (4.6), the 
n \ .5 ** i inutt" 1 <sre related through eqns. (2.19) to (2.2l). 

If T.er clear that for a given tube radius n^ and n^ 
h*v’ r of u f ixed for optimal cooling. Therefore, for a given 
tub 1 * 1 mdiuc oqn. (4.5) and (4.6) are coupled differential 
oq.j.-i n r > in i.vo variables, namely, n 2 and n 3 . 

4..' if! J i IhL SGhDITIOKS 

3 . not- wo hove only 2 variables, we need to establish 
jn.ti.il conditions only for these two. In the radiation zone, 
which con be represented by x < o, we have saturation condi- 
tion viz. n 2 = -^ . Therefore, for our purpose 

n G 

“"2 » 

n 3G = n 3 | *< o = n 3 | x < 0 

The values of n 3Q have been printed out from the computer 
program represented by the flow chart given in Fig. 3.5. 




Table 4.1 


Initial Conditions for Various Tube Radii 


S .No . 

Radius ,R 
cm 

.™ ■ 

n2QXlO (m ) 

-20 

n 30 x10 1 

1 

0.1 

8.8 

1.13 

2 

0.3 

12.7 

0.46 

3 

0.5 

13.3 

0.29 

4 

0.7 

13.5 

0.22 

5 

0.9 

13.6 

0.17 

6 

1.0 

13.7 

0.16 

7 

1.1 

13.7 

0.14 

8 

1.3 

13.8 

0.12 

9 

1.5 

13.8 

0.11 

10 

1.7 

13.8 

0.09 

11 

1.9 

13.8 

0.085 

12 

2.0 

13.8 

0.08 

13 

2.1 

13.8 

0.07 




t 



0e U / £ U ( £ ]2>A2>) j.0 A^isuap SS2)]UO!SU3UJIQ 


Fig. 43 Dimensionless density of level 3 in the post radiation zone 
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Table 4.1 gives the values of n 2 Q and n^Q for various tube 
radii, corresponding to the optimum operating conditions. 

By choosing an appropriate step size, the dimensionless 
number densities, i.e., n 2 /n 2 Q an< ^ n 3/ n 30 ^ ave ^ een obtained 
as a function of using R-K-G method for all tube radii, 
and have been summarised for tube radii 0.1 cm, 0.5 cm and 
2 cm m Figs. 4.2 and 4.3. 

This exercise enables us to obtain the number densities 
as a function of “<2 which will be used m the next section 
for evaluating the cooling (or heating) effect in the post- 
radiation zone. 

4.5 ENERGY REMOVAL RATE IN THE POST-RADIATION ZONE 

As explained earlier m Chapter 2, the cooling or heating 
effect is brought about by a change m the translational energy 
of the gaseous mixture. For a system of CO-CO 2 , the rate of 
change of translational energy per unit volume can be written 
as given m eq. (2.16). As has already been indicated, the 
term involving k 21 is negligible for the CO-CO 2 system. 

Hence, eqn. (2.16) can be written as ; 
dE 

HIT = n 3 k 34^ E 3“ E 4^ ~ n 2 k 23( E 3~ E 2^ + n 3 k 32^ E 3~ E 2^ ( 4 * 7 ) 

substituting proper expressions for , 


dE 

( H3-E 4 )-n 2 C23n 4 (E3-E2) + n 3 C 32 n x^ E 3'" E 2^ 


(4.8) 
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Here n 2 and are functions of -XT • The cooling rate is the 
negative of the rate of change of translational energy. Eqn, 
(4.8) gives the rate of change of translational energy per unit 
volume at a given value of XT . Integrating eqn. (4.8) over 
the entire volume of the post-radiation zone, we can get the 
total rate of increase m the translational energy m this 
zone. Denoting the negative of this quantity by E^ 2 , we have 


where 


°o dE, 

E R2 ~ ” £ dt 

A dx 

^max 

dE, 

= -uA f 

0 


= uAI n 

^max 

dE, 

1 = — f ( 

n 0 



‘max 


that n 2 and n^ have almost dropped to zero 


(4.9) 


(4.9a) 


has been chosen such 
values at r C= T max » 


Once n 2 and n 3 have been calculated as functions of , 
as represented m Figs. 4.2 and 4.3 for representative tube 
radii, the values of I can be calculated by Simpson’s rule 
using eqn. (4.8), together with the initial conditions given in 
Table 4.1. These values of I n are listed m Table 4.2. It is 
seen that I is positive, which means that there will be a temp, 
drop in the post-radiation zone. This is calculated in the 
next section. 
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4.6 TEMPERATURE DROP IN POST-RADIATION ZONE 

Temperature drop *v©£ in the post-radiation zone can 
now be easily calculated by a simple energy balance : 


C(J ’ C p ) co + (fCphoJ ^ e 2 = E R2- 


or 

^Qq 

z. 



[( J*c ) + 

k v p'co 


(f C ) 

p'cor 


(4.10) 


Using the values of given m Chapter 3, ^©2 has been 

calculated, and is tabulated in Table 4.2, as a function of 
tube radius at the optimum operating conditions. 


It is seen that there is a temp, drop of the order of 
2 to 3K m the post-radiation zone, which can be added 
to the temp, drop, A 9-^ m the radiation zone. It may be 
pointed out thatA ©2 is independent of the flow velocity, 
whereas is inversely proportional to the flow velocity 

(see Table 3.3). 
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Table 4.2 

Temperature Drop in the Post-Radiation Zone 


Tube 

Radius 

(cm) 

Integral I 
(J/m 3 ) 

a© 2 ( k ) 

0.1 

0.47 

2.0 

0.3 

0.49 

2.92 

0.5 

0.38 

2.54 

0.7 

0.36 

2.34 

0.9 

0.37 

2.45 

1.0 

0.38 

2.5 

1.1 

0.39 

2.57 

1.3 

0.40 

2.66 

1.5 

0.40 

2.74 

1.7 

0.42 

2.82 

1.9 

0.43 

2.89 

2.0 

0.43 

2.90 

2.1 

0.43 

2.90 
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Equation (3.1) 


can also be rewritten as 


n. 


- t!r- + 1^-1^21^23^ [^34 +k 32^ k 34^ “ k 32 k 


12 


^12^ 1 34 +k 32 +k 34^ 


23 


(5.2) 


ror 't,uura aon should be equal to n 2 . In the expression 
given above the first term is equal to one because W 12 is equa 
to ^21* Therefore, to make n^ equal to n^ , , the second 
expression on the right hand side of E^eqn. (5,2) should be 
made negligible as compared to 1. This is achieved by making 
the denominator sufficiently larger in comparison to the 
numerator. This condition is achieved when 


^12^ A 34 +k32+k 34^ ^ 10 ^^21 +k 21 +k 23' ^ A 34 +k 32 +k 34^ 10k 32 k 23 

(5.3) 


Ncq looting k 21 for reasons explained in Section 3.8 

substituting proper expressions for rate constants, 

condition for saturation is given as 

/ v .^n -i n A 

W 12 = i°(^i +c 23 n 4^ ~ aT7+C 


34 n "^32 n l +C 34 n T 


and 

the 


(5.4) 


The most conservative estimate for laser power calculation 
is the one when levels 2 and 3 are assumed to be empty. 
Physically this situation occurs only at the beginning. This 
fact further justifies dropping of the radiative trapping 
terms. At the beginning, since levels 2 and 3 are almost 
empty, F 2 ^ and are equal to 1, leaving A 21 and Ag 4 


unaltered 
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> <x ' h '-hi 


o sunolification, n 4 could be replaced by n. 


l H 


3 Tlj _| X S 


and nj by n«. Further, n G is equal to r times n T and 
equal co (]-r)rij. IJith these modifications, the saturati 

contn ti on becomes : 


on 


W i o 


-- 10(. ?1 ,-c 23 (l-r)n T -) 


10(C 32 C 93 ) ( r) ( l-r)n^ 

J ( rl ^32 + ^ ^ ) C 3 ^ ) O - j , 


(5.5) 


In eqn . (5.5) if r and n T are known, vV l2 can be deter- 
mined .-a mg the other constants that have been derived in 
Sections 2.4 and 2.5. The values of r and n^ are the operating 
onoc, vj'n : ch have been derived in Chapter 3 for various tube 
radii and summarised in Table 3.2. 

fk.o intensity of the laser beam required is given by 

the express 1 on [ 15] 


^ 12 [ ^ 2 ~‘^1 ] n J-_ 

K g G n i 


(5.6) 


where k is given by 


. > 2 &21 
n = 2i 
~ 8 % 


(5.7) 


If the { vu 1 tipli cities of levels 1 and 2 are assumed to be the 
same. For CO gas, corresponding to levels in Fig. 2.1, 

a 2143 err and A 01 = 33.8 sec" 1 which gives us 

>21 21 


-11 


2 -1 
m sec 


2.928x10 




(MU4 U I) JdMOcj 
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V 


which is the line shape function, is given by 



(5.8) 


for the case of Doppler broadening. 

has been explained in Section 2.7 and for CO gas, the 

Q 

value has been comouted as 1.507x10 Hz. 

All these values put together m eqn. (5.6) gives us the 
intensity I is multiplied with the area of cross section of 
the beam, which is the same as that of the tube through which 
the gas flows. The results of the calculations are shown in 
F ig . 5.1. 

5.3 COMPARISON WITH CONVENTIONAL SYSTEM 

In this section, the proposed system is compared with a 
conventional system operating at similar conditions of total 
pressure, partial pressures, flow speed and tube radius. We 
evaluate the temperature difference required across the tube 
wall, for the same amount of heat removal rate as for the 
proposed one. 

For this calculation, we assume that the convective film 
at the outer surface of the wall and tube wall, do not provide 
any resistance to heat flow. The resistance to heat flow 
therefore may be calculated on the basis of the convective 
film on the inner side of the tube only. 
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'Jg f li'U check whether the flow is laminar or bur- 
wulenc, to decide on the appropriate correlation to be 
used tor evaluating the heat trans Ter coefficient. For 
this purpose we evaluate the Reynolds number, given as : 



where 

u is the flow speed (m/s), 

D is the diameter of the tube in which the flow 

occurs (m) 

"“*1 

and p is the viscosity of the gas flowing (kg m S ) 


The viscosity is evaluated using the expression [ ] : 


M- 


2.669x10 


6 


o^-CL 1 1 


(5.9) 


where 


cr 

and 


is the molecular weight of the gas, 
is the absolute temperature (K) , 

is the characteristic diameter of the molecule (A ) , 
is a slowly varying function of the dimensionless 
temperature . Both and -A are tabulated 

functions . 


Equation (5.9) is found to be applicable to monoatomic as 
well as polyatomic gases [13 ]• 
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For CG at T = 300 °K 
cr --= 3.39 R 

^ ~ 1 • 064 
•F ri j ‘ 'T ■ = GO gives 

M- - 1 . 73 x 10 ""“* 
7 :r J 0 o it T - 300 °K 

o = 3.996 (R) 

• iKj Jl = 1.233 

M* 

- v. i ' for ,7 - 44 gives 
p = 1 . 5 xlO ~ 5 


:osi ty for a mixture of gases is evaluated using [ 13] 


n X i ^i 

i g x ~|.‘ 

i 0 iJ 


(5.10) 



n is the number of chemical species in the mixture, 
and Xj are the mole fractions of species i and j, 
p. and p are the viscosities of species i and j at 

i J 

the system temperature and pressure 
and and I-Ij are the corresponding molecular weights. 



O j 


O'.tV C >Zrj T t 


, I’ > same or 


he viscosity for CO and C0 2 are of 


f ‘ nrt ' TOle faction of 


C0 2 is very small, we 


or 'm j n 


"' iV - mu l • jj. by u „ . 

,:1X co ’ Using this approximation, and a 


1 ; i/s for 


a tube radius 1 cm 


t\ ^0.015 


*>'V UGn sity (J» ) has been evaluated using ; 


* ■ ^co * -f co2 


(5.12) 


-to dnd ^co2 aro calcu lated using eqn . (3>10) an(J for this 

particular case, their values are • 

^0 0 = 126,6 x 10-6 ^g/m 3 


11.6 x 10~ 6 kg/m 3 


The Reynolds number for other tuh» i n 

uoe raai1 will also be of the 

same order. For a flow throuqh a Hurt , ■> ^ „ 

y a auct » this value of Reynold 

number indicates that the flow is laminar. We can now use an 
appropriate relationship to find the heat transfer coefficient. 

For a fully developed laminar flow through a duct 


N = 4.364 
u 


(5.13) 


for a constant wall heat flux [l]. N(j> the Nusselts number is 
given by 





who 'G 


h 

k 

and D 

The thermal 


where 



M 

and \x 

For CO, 

M = 

c p = 

and 

P- = 


which gives 


and for CO2 


M = 




o 

is the heat transfer coefficient ( W/m K) , 

is the thernai conductivity of the gas (W/mK) 
is the diameter of the Lube (m) . 

conductivity k is calculated using [13], 

k = (C P + 1 i )l1 (5 - 14) 

is the specific heat of the gas (j/kg K) , 
is the universal gas constant and has a value of 
8314 (j/kg.mole K) , 
is the molecular weight, 

is the viscosity of the flowing gas (kg m ^ s . 


28, 

296 ‘ 9 (“kg“id 


1.78xlO~ 5 


0.026 


W_ 

mK 


48 


188. 9( 


J 

kg°K 


) 



and 


1.5X1Q- 5 jg 


which gives 


'co 


0.0161(W/mK) 


The thermal conductivities of gas mixtures at low Density may 
be estimated by a method analogous to that previously given 
for viscosity (eqn. (5.10)) [13] : 


n x . k . 

v — £ 

mix is=1 n 

S x 

j=l 3 13 


(5.15) 


The terms m the expression given above have the same meaning 
as for eqn. (5.10). 

Since the values of thermal conductivity k, for CO and CO 2 
are almost of the same magnitude, one could use, 


k 


mix 


0.024 


w_ 

mK 


Using this value of k, we can write 

hD = 0.1047 (5.16) 

using eqn. (5.13) 

The temperature difference £»T, required to be maintained 
across the tube wall for effecting a heat removal rate, 


can be calculated using ther -expression 



E R1 hS A A>T 


(5.17) 


where 

is the surface area across which heat flow 
2 

occurs (m ) . 

Surface area S A for a tube of diameter D and length L is 
ttDL. Using L D, as in the proposed laser induced cooling 
scheme, AT can be calculated for the various values of 
given in Table 3.2. For example, for R = 1 cm, we get 


and 


h ^ 5W/m 2 K 
T 10~ 2 K , 


Thus we see that the temperature differential required between 
the coolant and the fluid stream to be cooled is indeed very 
small for removing the quantities of heat that a laser induced 
cooling system removes. 

However, the value of heat transfer coefficient obtained 

through eqn. (5.13) appears to be unrealistic for this case. 

This needs to be improved upon. Perhaps eqn. (5.13) is not 

applicable for flow of gases through a duct at pressures of 

the order of 100 mtorr, as is the case here, especially when 

—2 

the duct diameters are also of the order of 10 m or so. 
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This comparison was nictde on the basis of heat removal 
rates. It should, however, be noted that in laser induced 
cooling system, one can achieve higher temperature drops over 
very small axial distances by reducing the flow speed. Such 
large temperature drops over a small distance cannot possibly 
be achieved m a conventional system. 



CHAPTER 6 


CONCLUSIONS AND SUGGESTIONS 

In this work, steady state cooling of a mixture of CO 

and CO2 , using laser radiation has been considered. The 

mixture flows In a tube and the laser beam falls perpendicular 

to It. The optimal operating conditions (partial pressures 

of the constituent gases) for maximum cooling have been obtained 

for several tube radii. These, together with the cooling rates 

are listed in Table 3.2. Although cooling will take place of 

if the total pressure is less than 200 rat on,- [ ], optimal 

operating pressures are m the range of 90 - 100 rntorr for tube 

radii ranging from .1 cm 'co 2 cm. The maximum volumetric 

3 3 

cooling rates achieved are m the range of 6 W/m to 90 W/m , 
and are higher for low tube radii. 

The Important effect that has been included is the radia- 
tion trapping. It is because of this phenomenon that cooling 
rates are high for low tube radii. With an increase in the tube 
radius, the cooling rates come down significantly. However, for 
larger radii, the volume Integrated cooling rates are larger. 

For a tube of diameter 1 cm, the volumetric cooling rate is 

O 

25.8 W/m and the volume Integrated cooling rate is 14.5 pW. 
These values are obtained for a total pressure of 92.7 rntorr 
and a partial pressure of CO 82.6 rntorr. A temperature drop of 



12° can be achieved ”/i tn a flow speed 0.1 m/s in a tube of 

diameter 1 cm. 

In this work, analysis has been done for flow in a tube 
with a circular cross section illuminated by a laser beam 
which is also circular m cross section. Analysis for 
rectangular and other sections of tube as well as laser 
beam could be done on similar lines. Further, the total 
cooling rate per unit laser power could be maximised and 
optimal operating conditions determined. The dependence of 
the rate constants on temperature could also be accounted for. 
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